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Abstract We discuss the history of X-ray observations of the ?/ Car system, con- 
centrating on the periodic variability discovered in the 1990s. We discuss the inter- 
pretation of these variations, concentrating on a model of the system as a ‘'colliding- 
wind” binary. This interpretation allows the physical and orbital parameters of 1 / Car 
and its companion star to be constrained. 


1 The 5.5 Year Period 


Periodic variability is a key diagnostic in astrophysics. Astronomers long sought 
periodicities in rj Car's behavior ever since its giant eruption. Before 1870 Rudolf 
Wolf and Elias Loomis suggested photometric cycles of 46 and 67 years, respec- 
tively 1 99, 60]; a century later Payne-Gaposchkin proposed a cycle of 15-16 years 
[73], and Feinstein & Marraco suspected a 3-year periodicity [35]. The complex 
behavior of the observed variations from rj Car in broad-band photometry and in 
spectrographic observations made determination of the true period difficult. 

The discovery of 17 Car’s first confirmed period depended on observations span- 
ning twenty years or more. This began with the recognition of a so-called “spec- 
troscopic event” by Zanella, Wolf and Stahl in 1981 [100, 95]. These “events” rep- 
resent a dramatic and complex change in t] Car's visible-band spectrum in which 
high-excitation He I and [Ne III] emission lines abruptly fade, then recover. Sim- 
ilar changes happened in 1948 and 1965, which suggested to them a recurrence 
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timescale of 17 years [100]. They were on the right track, but did not know that 
other such events had occurred (as we deduce today) in 1953, 1959, 1970, and 1975. 
Near-infrared photometry from South Africa obtained by Patricia Whitelock and her 
collaborators[96] over 156 nights in 1975- 1994 showed a significant 5-6 year cycle, 
but this cycle was confused with other near-infrared variations on other timescales. 
Meanwhile new spectroscopic events occurred in 1986 and in 1 992, though the for- 
mer was not widely reported. 

The breakthrough was provided by a set of 27 measurements of the He I 10830A 
emission line obtained by Augusto Damineli at Brazil’s National Astrophysical Lab- 
oratory'. These observations, augmented by 4 additional measures from other ob- 
servers from 1981-1 995 , showed a clear, rapid weakening of this He I line every 5 .5 
years, lasting only for a few weeks (a duty cycle of < 5%). Damineli showed! 16] 
the fading of the He I 10830 line occurred at nearly the same time as brightening 
of the system in the near 1R, and that the 5.5 year period recovered not only the 
“spectroscopic event” seen by Zanella, Wolf and Stahl in 1981, but other "spectro- 
scopic events” observed in the previous 50 years. Damineli successfully predicted 
the occurrence of the next minimum at the end of 1997. 

Thanks to multi -wavelength observations of i] Car during the subsequent minima 
in December 1997 and June 2003, we know that the fading occurs throughout nearly 
the entire electromagnetic spectrum from the radio through the 1R, optical, UV and 
X-ray region (the situation in the Gamma-ray band is less clear). The period is al- 
most exactly the same in every wavelength region[21]. Detailed monitoring obser- 
vations since 1997 have refined the period to 2024 days (5.54yr) with an uncertainty 
that appears to be no more than a week and may be as small as ±2 days[20]. 

What regulates the cycle? Following Zanella et al., Damineli initially supposed 
that each spectroscopic event is a single-star “shell ejection”, a temporary large 
increase in the stellar mass-loss rate, reminiscent of phenomena seen in S Dor vari- 
ables and Luminous Blue Variable stars (see, e.g., [46)). But ij Car’s variability cy- 
cle is too long to represent a pulsation or rotation period, since those periods should 
be only days or weeks. A thermal/rotational recovery cycle could have a timescale 
of several years [31], but this seems unlikely to produce so regular a period. The 
He I 10830 variation seems strictly periodic and strongly suggests a simple clock- 
like explanation: a companion star in orbit around rj Car. If the 5.54 year period is 
assumed to be the orbital period of the system, the high total mass (probably 150 
to 200 M. ) implies a semi-major axis close to 17 AU. The short duty cycle implies 
that the phenomenon occurs near periastron passage in a highly eccentric orbit. Pro- 
posed orbits [17, 18, 23, 24, 25, 26, 53] are mostly comparable in size and shape to 
the orbit of Halley’s Comet. 

But the UV-to-IR emission regions include both the wind of rj Car and the de- 
tritus of the eruptions from the 19 1 * 1 century and are spatially complex. Emission 
from one region may be absorbed by cold gas in another. Dense clouds may shadow- 
one parcel of gas from the photospheric radiation but not its neighbor. More dis- 
tant gas reacts more slowly to changes iti the inner system because of light-travel 
time, which may be days or even weeks. These complications make interpreting the 
spectrum at these wavelengths difficult. 
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A binary system composed of luminous stars must produce significant amounts 
of X-ray emission resulting from the high speed collision of the massive stellar 
winds[5,76| which constantly flow off the stellar surfaces. This wind-wind collision 
will produce a “bow shock” where the wind kinetic energy is converted to thermal 
energy which radiates mostly in the X-ray band (as discussed below). Such emission 
is an important probe of massive binaries since it originates in a localized region 
(defined by system parameters like the separation between the two stars and the 
relative mass loss rates and wind speeds), and it penetrates through large columns 
of intervening cold material. What do X-ray observations tell us about rj Car? 


2 X-ray Emission and Variability 


Studies of X-ray emission from rj Car go back 30 years or more, nearly to the begin- 
ning of X-ray astronomy. The earliest data were obtained by non-imaging detectors 
with fields of view defined by collimators, flown on rockets or orbiting platforms 
[42, 39, 36, 77, 69], They established i; Car as a source of X-rays at energies E > 5 
keV, and detected iron K-line emission at 6.7 keV. They also provided hints of vari- 
ability, though direct comparisons of X-ray flux were hampered by differences in 
instrumentation (fields of view, sensitivities, etc). 

Today we know that many X-ray sources, both discrete and diffuse, exist near 17 
Car, so it was not until the advent of imaging X-ray telescopes that X-ray emission 
from 7] Car could be cleanly disentangled from this contamination. The first X-ray 
images of rj Car and the Carina Nebula were obtained in the early 1980s with the 
HEAO~2/Einstein satellite observatory [78, 79, 6[. Those observations showed that 
at energies above 3 keV, )/ Car was the only source visible in the field. They showed 
that X-ray emission from r\ Car consists of at least two components: an extended 
(~ 20") “soft” ( E < 1 keV) component (likely produced by the collision of material 
from the 1840’s eruption with interstellar material) surrounding a point-like, “hard” 
(E > 3 keV) heavily absorbed component centered on the optical position of rj Car. 
This “point-like” component has a temperature of 4 — 8 x 10 7 K and its X-ray emis- 
sion is absorbed by intervening material having a column density of 5 — 50 x 1() 22 
atoms cm 2 [6. 12, 40]. The observed luminosity of the hard, point-like source mea- 
sured by Einstein (adjusted to the current best estimated distance of 2300 pc) was 
ss 4 x 10 33 ergs s~' in the 0.2-4 keV band; correcting for absorption, the luminosity 
of the source is about 2 x 10 34 erg s -1 . While these observations were key to estab- 
lishing the overall character of the high-energy X-ray emission from 17 Car, these 
relatively short “snapshots” could not really address issues of temporal variability. 

Figure 1 shows a Hubble Space Telescope image of 17 Car and the Homunculus 
nebula superimposed on a more recent X-ray image from the Chandra X-ray Obser- 
vatory. This image clearly shows the extended “shell” of X-ray emission first seen 
by Einstein surrounding the Homunculus along with the point-like X-ray emission 
coincident with the optical position of the star. The point-like, hard emission was 
initially interpreted[79, 6[ as the collision of a fast wind from i| Car with a circum- 
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stellar “shell” about 1000 AU from the star. As we shall see, the temperature of the 
emission requires a wind velocity of 2000-3000 km s" 1 or more if it's solely pro- 
duced by conversion of the wind kinetic energy to heat. This is about 4 times faster 
than i/ Car’s bulk wind seen in visual-wavelength spectra - a highly signficant point 
in retrospect, though it was not emphasized at the time of the Einstein observations. 
The X-ray absorbing column depth was thought to be consistent with absorption in 
a 3 solar mass cold shell, roughly agreeing with what was known of the mass of the 
Homunculus at that time. 


2.1 Discovery of X-ray Variability 


X-ray variability of rj Car was first detected by the ROSAT X-ray observatory in 
short observations obtained in 1992-1993 analyzed by G. L. Rawley. She showed 
that the central source was not visible in June and August of 1 992 at energies above 
2 keV, but was clearly detected in January 1993 [7]. The flux above 2 keV varied by 
a factor of two in a 4-month interval , requiring a tripling of the mass-loss rate if the 
X-rays were produced by collision of the wind with a circumstellar shell. However, 
variations in the temperature of the emission or changes in the intervening absorbing 
medium could not be ruled out. 

Luck often plays a role in astronomical discovery; the mid-1992 ROSAT observa- 
tions just happened to coincide with one of the spectroscopic events, and to coincide 
with the first observations of the star|33] by the radio telescopes of the Australia 
Telescope National Facility. At the same time t; Car faded in the hard X-ray band, it 
also faded in the radio [33], but brightened in the near-IR [96, 16). Thus, when the 
5.5-year period was identified by Damineli shortly thereafter, an obvious question 
was whether X-ray variability was another aspect of the same cycle. 


2.2 The X-ray Variability Cycle 


Monitoring observations with the Proportional Counter Array (PCA) on the Rossi Ar- 
ray Timing Explorer (RXT£)[4| provided the next X-ray breakthrough. Unlike other 
X-ray observatories, RXTE can view q Car at any time of the year making it a near- 
ideal instrument to study the variability of the 2-10 keV flux. RXTE observations of 
rj Car began in February 1996 and have continued through 2009(52, 14]. Ty pically 
RXTE obtained 2-4 observations of ij Car per month, with daily observations during 
important intervals. The RXTE X-ray count rates vs. time are shown in figure 2. 

The RXTE observations show the important temporal signatures of the 2024- 
day X-ray variability cycle. For about four years out of the 5.5 year cycle, the 2- 
10 keV X-ray luminosity is in a “quiescent” state, relatively constant at a level of 
~ 4 x 10 34 ergs s~‘ , but with erratic small-amplitude variations. About a year be- 
fore the fading of the He I 10830 line, both the intensity and the irregular variability 
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I’ig. 1 t; Car in X-rays. The upper figure shows an X-ray image from the Chandra X-ray ob- 
servatory (in green) combined with an optical image from the Hubble Space Telescope (in red). 
The unresolved white point source at the center of the image, at the position of the optical star, is 
the source of hard X-ray emissions. The lower figure shows five Chandra AC1S-I images of 1 / Car, 
highlighting the variation of the central pointlike source. The streaks at various angles to the central 
source are detector artifacts. 


progressively grow, and we see conspicuous X-ray “'flares” vvi th durations of weeks 
or days, typically two or three months apart. At the end of this phase, the average 
X-ray luminosity in the 2-10 kcV band has roughly doubled to reach a maximum 
luminosity of « 10 35 ergs s" 1 , while the interval between flare peaks drops to less 
than l month. About 40 days after reaching its maximum, the emission swiftly (with 
a few weeks) falls to a level below KXTE's detectability. This rapid decline occurs 
along with similarly rapid fadings of the He I 10830 line strength[2i| and the IR 
brighlness[97|. The X-ray emission usually remains in this “low state” for about 3 
months. At this time the X-ray intensity is lower than or near the cosmic X-ray back- 
ground in the ~ 1° PCA field of view, though r/ Car is still detectable throughout 
the minimum by other X-ray observatories with sufficient spatial sensitivity! 1 3, 40| 
at luminosities ranging from ~ 6 x 10 32 ergs s~ 1 to 5 x 10 33 ergs s _1 near the end 
of the minimum|401. At the end of the low state the X-ray intensity rises rapidly to 
the “quiescent" level. 

X-ray spectra of i) Car show that the observable emission associated with the 
star is mostly in the 2-10 keV band, corresponding to a temperature range T = E/k 
of roughly 20—1 10 million K. The PCA spectra show little flux above 10 keV at any 
phase. Near X-ray minimum, the observed emission declines more at lower energies 
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0.0 0.5 1.0 1.5 2.0 



Fig. 2 The RXTE X-ray “lightcurve” of )j Car at energies above 2 keV and below 10 keV. as 
measured by the PCA on RXTE. Three X-ray minima can be clearly seen near phases 0.0, 1 .0 and 
2.0. The “flat” bottom of the X-ray minimum seen by RXTE is produced when the X-ray flux from 
i) Car falls below the cosmic background level (shown by the horizontal dashed line). The blue 
and red lines show the RXTE measures offset by 2024 and 4048 days, respeclively.Adapted from 
Corcoran etal. [15], 


showing the build up of circurastellar absorption, and the spectrum remains quite 
hard for an extended period even after the minimum ends. 


3 Colliding Wind Emission as a Diagnostic of Massive Binary 
Systems 


If the pan-chromatic periodicity is produced by the orbit of a companion star around 
>1 Car (as seems likely) then this offers the unique opportunity to dynamically mea- 
sure the mass of a star near the Eddington Limit using conventional binary-star 
radial velocity analysis. Unfortunately attempts to solve the radial velocity curves 
of observable emission lines in the UV-to-IR [17, 23, 19, 27, 68j have not given 
trustworthy orbital elements, mostly because these lines are formed in a spatially 
complex mass flow and consequently are broad, asymmetric, variable, and highly 
model-dependent. For this reason, the X-ray variability, produced by the orbit-driven 
variation in the wind-wind collision zone, provides a very useful (though indirect) 
diagnostic of the orbit. We briefly review the physics of X-ray generation in collid- 
ing wind binaries to see why this is so. 

Colliding wind binaries (CWBs) produce thermal X-rays in the strong shock 
where the stellar winds collide. The relevant geometry' is crudely sketched in Figure 
3. The shock is concave toward the star with the “weaker” wind. i.e. the star with 
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the lower wind momentum flux. The vertex of this shock, where the winds collide 
head-on close to the line of centers, is the “stagnation point” where the winds are 
fully decelerated. In the idealized case a contact discontinuity separates the shocked 
stellar winds from the primary star on one side and secondary star on the other. 
The temperature of the shocked gas on each side of the contact discontinuity is to 
first order determined by the pre-shock velocity of the stellar wind from the star on 
that side. This is important for q Car since the wind speed of q Car itself is only 
~ 500 km s 1 , which would produce only soft X-rays which would be unobservable 
due to heavy absorption by the cold gas in the Homunculus. The fact that hard X- 
rays are observed from q Car indicates the presence of a fast flow in the system. 
The main X-ray emission processes are thermal Bremsstrahlung from hot electrons, 
and inner-shell (K, L, M) transitions in nearly stripped ions. Standard examples of 
CWBs are HD 193793/WR 140, y 2 Velorum, and HD 5980 in the SMC. These 
massive binaries have observed X-ray luminosities in the range 1 0 32 — 1 0 36 ergs s'~ ! 
« 0.02-300 Lq, varying on orbital timescales. Though large, these luminosities are 
of course only a millionth or so of the immense stellar bolometric luminosity. 



d2 


Fig. 3 A sketch of the q Car colliding wind binary system. Winds from two stars flowing radially 
from the stellar surface collide and form shocks separated by a contact discontinuity. In this simple 
model, the shocked gas fonns two curved shock fronts between which the heated plasma flows 
out along the contact discontinuity. Temperatures on either side of the contact discontinuity are 
determined by the velocities of the stellar winds just before the collision. 

Theoretical descriptions of X-ray emission from wind-wind collisions in mas- 
sive binaries go back to the 1970’s, and the details can be quite complex [5, 76]. 
The essential behavior, however, can be roughly described by the following order- 
of-magnitude analysis by K. Davidson (private communication). Shocked gas cools 
by expansion and by emission of radiation as it flows outward away from the stag- 
nation point. Characteristic flow speeds in that region should be comparable to the 
isothermal sound speed w = (P/p) li ' 2 oc T i/2 . Based on Fig. 3, we also expect 
the lateral radius of the hottest region to be similar to A. the distance between the 
secondary star and the stagnation point. Therefore the expansion cooling time is 

^ d-2 ^ do 

4xp ~ ~w ~ (10 4 cm s~'~K C 2 ) F'/ 2 ' 
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Meanwhile the radiative cooling time is approximately the thermal energy per unit 
volume divided by the Bremsstrahlung emission rate per unit volume: 


J>n e kT _ CT [ ' 12 
AT'/ 2 n| n e 


( 2 ) 


where n e is electron density, A » 10" 27 ergcm 3 s 1 K '/ 2 . and C. w 10 1L5 cm 3 sK l/2 . 
Thus 

^ « (10-“crK)^. (3) 

trad ' 

We can relate T and n e in the hot shocked gas to the speed and mass loss rate in 
the secondary wind. According to the Rankine-Hugoniot adiabatic jump conditions 
across a shock [58, 80|, the maximum speed of sound and maximum temperature at 
the shock front are 

wWx = (3/ 16) 1,/2 v and T mm = w 2 ax , (4) 


where v is the pre-shock wind velocity, mu is the proton mass and a is between 0.6 
and 1 .0 for the expected range of chemical composi tion. If the relevant average T in 
the above equations is half the maximum temperature, we find 


T « (10 7 K) 


1000 km s 


(5) 


As noted above, this is the reason why q Car’s primary wind does not pro- 
duce observable X-rays: its typical speeds (which, it turns out, are latitudinally 
dependent^ 1J) of 200-1000 km s 1 lead to average temperatures of only a few 
million K, and the resulting soft X-rays (E < 0.7 keV) are absorbed by the large 
column of gas {N H ~ 5 x 10 22 cm ~ 2 ) in the Homunculus nebula. The observed 2- 
10 keV X-ray spectrum indicates T > 4 x 10 7 K, which requires v > 2000 km s 1 . 
So far as observable X-rays are concerned, q Car’s dense primary wind merely pro- 
vides a wall for the secondary wind to collide into. 

The jump conditions across the shock require that n e in the shocked gas is about 
4 times as large as the electron density of the pre-shock wind . The pre-shock density 
is related to the secondary star’s mass-loss rate by mass conservation: M = 4nr 2 p v. 
If the distance r from the secondary star to the shock is roughly equal to cA, and the 
pre-shock velocity of the secondary's wind is vz, then the density in the heated gas 

(It) 23 ' 2 g ' 1 ) M 2 

n e ~ -5 • ( 6 ) 

dz V’2 

Substituting for n e and vi in the earlier expression for r e xpArad, we find 


texp 

trad 


I0 / 


Mj 


Mr. yr 


1 AU 

d~> 


10 6 K 


3/2 


T 


(7) 
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Using M% ~ 10 ' >r> V/ . yr 1 |74], d 2 ~ 1 AU near periastron, and T ~ 5xl0 7 K 
(consistent with the observed 2-10 keV continuum), yields texpArad ~ 0.1. This 
means that expansion dominates the total cooling rate but not by an enormous factor. 

More formal and detailed analyses of colliding winds have been calculated [86. 
92,74), but the essential physics is captured by the simple analysis above. Numerical 
gas-dynamic simulations confirm our simplified reasoning in an order-of-magnitude 
sense, though incipient instabilities cause the shocked region to become spatially 
complex [74] . These detailed analyses suggest that the observed X-ray spectrum of 
) ] Car requires v 2 ~ 3000 km s _1 because of rapid cooling of the hot gas [74|. 

If fexp /trad A 1 radiative cooling dominates, and the shock structure becomes 
unstable [86[. For rj Car this is expected to occur near periastron in the shocked 
primary-wind [74[ on the left side of Fig. 3. In some proposed scenarios it occurs 
on both sides of the contact discontinuity, so that both shocks collapse near perias- 
tron [30, 61 1. In this connection it is worth noting that here we have ignored likely 
inhomogeneities, which may exist in the wind or may be caused by shock insta- 
bilities. If the density is inhomogeneous, then the meaning of n e in the expression 
for r e xp Arad becomes less obvious. On average, inhomogeneities tend to increase the 
fraction of input energy that is radiated. 

If the ratio £ — fisxpArad is small, then it approximately represents the efficiency 
for conversion of the secondary wind's kinetic energy to radiation. We noted above 
that the observable X-rays (E > 1 keV) in a CWB model for ij Car are powered by 
the companion star’s fast wind, whose energy outflow rate is A A v|/2. Consequently 
the observed X-ray luminosity should obey these approximate proportionalities: 


L x « £ x Kl 2 v 2 2 


M 2 

dTf^ 


X M 2 V 2 2 oc 



( 8 ) 


assuming that a constant fraction of the wind encounters the shock. 

A key point is the inverse dependence of L x on d 2 ; if the companion star fol- 
lows an eccentric orbit, then the varying distance modulates the X-ray luminosity 
such that the intrinsic X-ray luminosity emitted by the shock is highest near peri- 
astron and lowest at apastron. Therefore the modulated intrinsic X-ray luminosity 
can in principle be used reconstruct the orbit by determining the stellar separation 
as a function of time. However it is often non-trivial to derive the intrinsic X-ray 
luminosity from the observed luminosity, which includes the effects of variable ab- 
sorption from one (or both) of the highly structured stellar winds. If we exclude 
the complicated interval near periastron, the apparent flux variation shown in Fig. 
2 is dominated by the combination of the \jd 2 variation of the intrinsic flux and 
the phase-dependent absorption along our line of sight. These two effects depend, 
in turn, on three classical binary-star parameters: (l) the inclination i of the or- 
bit plane, usually assumed close to 45°, so that the orbital angular momentum is 
co-linear with the bipolar axis of symmetry-' of the Homunculus nebula [29[; (2) ihe 
orbit eccentricity e, which determines r(t)\ (3) the azimuthal orientation of the orbit, 
specified by the longitude of periastron to. customarily defined so that at periastron 
the secondary star is moving away from us at to = 180° and is on the far side of 
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the primary at co = 270°. As noted above, the intrinsic X-ray temperature can be 
used to constrain the speed of one (or both) of the stellar winds. The intrinsic X-ray 
temperature should vary little, because it depends chiefly on T and therefore on the 
wind speed v 2 , which we expect to be more or less constant. 


4 Modeling the RXTE X-ray Lightcurve of rj Car 

The first attempts [51, 52, 88 1 to explain the RXTE X-ray lightcurve of >j Car with 
colliding wind models were reasonably encouraging. Simple analytic models could 
describe the overall variation in X-ray flux and in Nh and the observed constancy 
of temperature T with plausible assumptions about the stellar winds, component 
masses and orbital inclination. The abrupt change from X-ray maximum to mini- 
mum is suggestive of an eclipse of the X-ray source region by the primary- star’s 
wind, so these first attempts assumed that the X-ray minimum was due entirely to 
an increase in absorbing column Nh, using simple, 1 -dimensional, spherically sym- 
metric wind models. Since the duty cycle of the X-ray “eclipse” is so small, it is 
very likely that the X-ray minimum occurs near periastron, when the stellar separa- 
tions and orbital geometry are changing most rapidly and column densities should 
naturally be large as the X-ray source becomes surrounded by the thick inner pri- 
mary wind. These initial attempts required a high eccentricity (e « 0.9. similar to 
the eccentricity of WR 140) and were consistent with an orientation in which the 
companion star is in front of the primary near apastron and behind at periastron. 
However, because the orbital geometry changes so rapidly around periastron, the 
duration of the X-ray minimum in these simple “eclipse” models was always far 
shorter than the observed duration. There are many modifications possible to cor- 
rect this, among them: 

• The parameters i and co can be altered [26, 30, 53, 70 J; 

• Tidal and other effects may suddenly increase the primary' wind’s density near 
periastron [23,26,81. 11,61]; 

• The large-scale structure of the primary wind is not spherical [8 1 1; 

• The orientation of the shock structure may be significantly distorted by orbital 
motion [74,70|; 

• Radiative cooling may destabilize the shock near periastron, possibly causing a 
general disruption or collapse of the X-ray source region |30. 61 ); 

• and other effects [87, 37,82,83], including radiative inhibition of the winds near 
periastron[72|. 

Most of these possibilities are not mutually exclusive; some of them appear quite 
likely , several may occur together, and some can reinforce each other. 

Throughout the cycle the observed X-ray flux depends critically on absorption 
in the intervening part of the dense primary wind; Nh obviously varies with phase 
in the orbit. The wind geometry is complex, since the secondary wind produces a 
cavity on one side of the primary w ind and orbital motion affects the cavity’s orien- 
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Fig, 4 Modeling the RXTE X-ray lightcurve of »j Car. Top: The solid blue line shows an 
estimate[52| of the X-ray light curve using an analytical model[92] with published orbital 
parameters! 17 1, assuming that the X-rays originate from a point near the apex of the shock cone. 
The dotted red line shows the results of the same analytic model based on more recently published 
orbital elements[70). Middle: Variations in column depth Np to the X-ray source along with the 
model variation|52|. Line colors are as above. Bottom: X-ray lightcurve (black) and model (red) 
based on a 3-D smoothed particle hydrodynamics mode![70] which better accounts for the varia- 
tion in column depth from the wind of the primary and successfully recovers the depth and duration 
of the X-ray minimum. None of the models account for observed cycle-to-cycle changes. The two 
upper panels are adapted from Lshibashi et a), j 52 1 .mid the lower panel is from Okazaki et al. |70J. 
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tation and shape. Analysis of the X-ray spectrum near the mid-point of the cycle|74| 
(i.e. nearapastron in the orbit) showed that the velocity of the companion's wind is 

« 3000 km s 1 , while its wind mass-loss rate is ps ICT 5 M 0 yr 1 . Recent efforts to 

model ?/ Car’s colliding winds with gas-dynamic simulations 1 74, 70, 72J attempt 
to account for variations in the large-scale wind geometry. Dynamical wind mod- 
els (711 and smoothed-particle simulations [70 j show the distortion of the wind 
by the secondary star’s orbital motion. Fig. 5, for instance, shows a smoothed- 
particle model of the wind-wind interaction at a number of orbital phases [70|. 
A post-periastron buildup of material on the forward-facing shock surface causes 
extra absorption that extends the eclipse. Thus a pure eclipse model with suitable 
parameters can account for the duration of the X-ray minimum (74, 70], but ap- 
parently only if the X-ray emitting region is small. If the X-ray emitting region 
is extended, the eclipse model has difficulty reproducing the extended duration of 
the X-rav minimum(72|. This does not disprove the complications and alternative 
models listed above, however; observations in other wavelength regions, beyond the 
scope of this review, are harder to explain than the X-ray minimum (61 ]. in sum- 
mary, at phases of the cycle not near periastron, i.e., not during a spectroscopic 
event, the models cited above are probably realistic in a semi-quantitative sense. 
They may be valid near periastron as well, but the complexity of the situation there 
makes this difficult to establish. 

Altogether, then, one can attempt to lit most of the X-ray “light curve” with just 
two orbital parameters, e and to. and the mass-loss rate and terminal velocity of the 
companion star’s wind. In 2000 Ishibashi estimated e « 0.8 and to « 200° by this 
method |53(. Implicitly using the same dependences, more elaborate calculations 
by Okazaki have led to e ~ 0.9 and to « 243° [70]. Thus, if the reasoning is valid, 
periastron occurs 20 3 to 70° before “superior conjunction” (the point where the 
secondary star is most nearly on the far side of the primary). Orientations in which 
the companion star is in front of q Car at periastron have also been suggested [34|, 
but these require either a very large buildup of absorbing material behind the X-ray 
source or a dramatic change in the shock emissivity near periastron passage. 


5 Dynamics of the X-ray Emitting Region 

Thermal colliding-wind emission is rich in X-ray lines that provide information on 
flow dynamics (via line widths and velocities), abundances (from line strengths) 
and temperatures and densities (from line ratios). The strongest lines are K-shell 
transitions in hydrogen-like and helium-like ions. In the case of rj Car’s wind, the 
observable X-ray features are broad lines of Si XIII, Si XIV, S XV, and S XVI, and 
Fe XXV at higher energies |41 1. The implied elemental abundance ratios appear to 
be near solar[8, 41 1, which is not surprising since neither rj Car nor its companion 
star is thought to be evolved enough to show enhanced silicon, sulfur or iron abun- 
dances at the surface. Unfortunately the carbon, nitrogen, and oxygen abundances. 
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which help trace the evolutionary state of the system, are all at energies below 1 keV 
and so cannot be studied this way because of circumstellar absorption. 

Observations in 2000-2003 with the transmission gratings on the CHANDRA X- 
ray telescope have resolved the broad lines and show that they apparently vary with 
phase in the 5.5-year cycle [41 1. Near times of X-ray minimum these lines shift to 
increasingly negative Doppler velocities while their profiles become broader and 
more distorted [3, 41 1. Interpretations of these changes are controversial. The ob- 
served variations may represent changes in the spatial distribution of the hot shocked 
gas [41 1. or perhaps the appearance of an additional outflow in the system [3 1 . So far 
only a single cycle has been observed so it's not clear if these changes are secular or 
phase-dependent. 

'The spectrum also shows a 6.4 keV K-sheil line of less ionized species of iron 
|8, 12, 40[. This is presumably fluorescence produced as the main X-ray emission 
passes through “cold” iron rich material, either in the wind of >] Car or in the Ho- 
munculus. The fluorescence line is broad, suggesting that at least part of it originates 
in fluoresced material in rj Car’s wind. During the minimum the line increases in 
equivalent width, but the measured increase is rather modest, which may mean ei- 
ther that the scattering zone is compact enough to be heavily obscured at that time, 
or else that the intrinsic X-ray emission actually declines during minimum. 

Thomson-reflected X-rays from the Homunculus have been resolved in CHAN- 
DRA images during the minimum when the direct X-rays are weak 1 13 1. Such ob- 
servations are important since X-ray reflection could (in principle) provide us with 
a 3-D view of the cold absorbing wind material as a function of orbital phase. Light 
travel lime delays are a complication, since light travel time for the reflected X-rays 
may be as much as 80 days towards the back-projected northwest lobe. 


6 The X-ray Spectrum During the Minimum 


The best measures of the X-ray spectrum during the 2003 minimum used the XMM- 
Newton and CHANDRA X-ray observatories [ 13.40,41 1. XMM's imaging detectors 
provide sufficient spatial resolution to exclude sources unrelated to >) Car, while 
spectral imaging with CHANDRA resolved the source even from contamination 
by the Homunculus and outer ejecta. These observations [40[ confirmed that the 
apparent X-ray flux decreased even at 10 keV. requiring either a column density 
Nh > 1 0 25 cm 2 lasting for two months or more or else a real decrease in the intrin- 
sic hard X-ray emission. Meanwhile, analyses of the A MM spectra |40] suggested 
that the relevant average column density increased by only a factor of 5 to 10 dur- 
ing the minimum, reaching a maximum of roughly 4 x 10 23 cm 2 . This was not 
sufficient to explain the observed hard X-ray decrease. Either the analysis under- 
estimated the maximum Nh, or else the intrinsic X-ray production considerably 
decreased, i.e. the event was not merely an eclipse [100, 25, 26] . 

During the minimum the hard X-ray flux varied strongly in timescales of days 
but not hours [40|; this seems reasonable if the lateral scale size of the emitting 
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region was of the order of 1 AU at that time. The minimum apparent flux in the 2- 
10 keV band was only 0.7% of the maximum seen by RXTE. The slope of the X-ray 
continuum above 5 keV does not vary much, indicating that the electron temperature 
of the hottest plasma does not vary significantly at any phase. The X-ray emission 
seemed to show two states [40] during the minimum: first a “deep” stale at the 
start of the minimum, in which the entire 2-10 keV flux was reduced, and then a 
“shallow” state beginning about midway through the minimum, when the 6 — 10 
keV emission increased and then remained constant for a brief period while the 
E < 6 keV emission continued to increase as the column density decreases. 

The 2003 minimum also revealed a previously hidden "central constant emis- 
sion” component of X-rays [40]. This dominated the system brightness at en- 
ergies E < 2 keV during the minimum, with an intrinsic X-ray luminosity of 
~ 10 34 ergs s _l in the 0.1-10 keV band. Its origin is unknown; if thermal, the 
observed temperature implies a flow velocity of 1000-2000 km s ~ ! . It may rep- 
resent the interaction of tj Car’s polar wind with nearby structures, perhaps the 
“Little Homuneulus”[54] (though it’s not clear that the Little Homunculus would 
survive any such collision for long); or it might be emission produced by outflow- 
ing shocked plasma interacting with cooler eircumstellar material. It might also be 
intrinsic emission produced by shocks distributed in one of the stellar winds. 


7 The Companion Star 


To date, there has been no generally accepted direct detection of the secondary star 1 * * * * , 
but for reasons noted below this is not surprising. In the absence of direct UV-to- 
IR information, the X-ray data provide perhaps the most useful (though indirect) 
constraints on its physical parameters. 

First, even if we were unaware of the minima and the 5.5-year cycle, the X-rays 
strongly suggest that jj Car has a companion star. The average X-ray to total lumi- 
nosity ratio, Lx/Lbol ~5x 10 6 , is very high for an evolved massive star, probably- 
indicating that some sort of companion object is needed to help generate the X-rays. 
There’s significant flux at high photon energies E > 5 keV, and such hard X-ray- 
emission in evolved massive stars is almost unknown except in magnetized objects 
(whose X-rays vary on short rotational time scales) and binaries. 

In a colliding-wind model, the X-ray temperature indicates the pre-shock speed 
of the wind; )/ Car’s X-ray temperature requires a wind speed of about 3000 km s“ l 
1 26, 74|. This is far higher than speeds observed in the primary wind. Indeed it is 
near the upper limit of wind speeds observed for any massive stars. For a radiatively 
driven wind from a hot (J^ > 21000K) massive star the wind terminal velocity is 


1 Iping et al. asserted that certain variable emission features in the 905-1 180 A wavelength range 

were directly associated with the secondary wind [48]. This may be true but it is controversial [44. 

61.81 1. The reported features were much narrower than one would expect for the fast secondary- 

wind, and in some models they can originate in the primary wind. The wind velocities derived by 

Iping et al . were also about a factor of 2-3 too low to produce the colliding wind X-ray emission. 
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related to the photospheric escape speed V esc as V x «3x V esc 157], so that a terminal 
speed of 3000 km s" 1 corresponds to a photospheric escape velocity of roughly 

1000 km s 1 . For a massive, luminous star the escape velocity is modified from 

the usual relation to acco unt for acceleration due to radiation pressure on electrons: 
Vgsc = l -T)/R, where M and R are the mass and radius of the star, and V 

the ratio of radiative Thomson acceleration to gravitational acceleralion(57]. 

In order to produce the observed X-ray luminosity, the same models show that 
the wind mass loss rate of the secondary star must be Ms ~ 10 0 M© yr 1 . This 
is exceptionally high for an O-type star, but seems marginally allowable for an ob- 
ject above 40 M©. The determination of M 2 is rather complex, depending on the 
poorly-constrained cooling ratio r e xp/f ra d discussed in Section 3 above. For exam- 
ple, inhomogeneities in the wind might increase t aKp /t n a- It’s possible that improved 
models can attain the observed X-ray luminosity with a smaller value of M 2 . 

In a ‘least surprising” evolutionary picture, one expects the companion ob- 
ject to be an early- or mid-O type star. Being less massive than j; Car it should 
be less evolved and therefore near the main sequence. A mid-O spectral type 
with 7j,jr ~ 35000 K would account for the photoionization seen in the ejecta 
[25, 26, 28, 43, 48, 47], Thus the most straightforward and most probable type 
of companion is a fairly “normal” star in the 30-60 M© mass range, though other 
possibilities can be devised. Other indirect arguments based on photoionization of 
circumstellar clouds by the UV flux from the companion star also point to a simi- 
lar, though perhaps somewhat more evolved, object [93, 62 1. The recent analysis by 
Mehner et al. [62 1 demonstrated that an 04—06 giant with L ~ 4 x 10 5 £,©, Teff ~ 
39.000 K, and M/ams 40-50 M© would fit the observed range of parameters quite 
well. Such a star radiates mainly in the hard-to-observe far UV, should be less than 
3% as bright as the primary at near-UV to 1R wavelengths, and produces no spectral 
features that would be readily detected when superimposed on the complex primary 
wind spectrum [44] . The lack of a direct detection of the companion is therefore 
unsurprising. The only disquieting detail is the companion's high mass-loss rate 
required in some X-ray models, which might be viewed as evidence that the com- 
panion is a Wolf-Rayet star instead of an unevolved C) star - making the system’s 
evolutionary status problematic. However, we noted above that the deduced value 
of Mi is fairly uncertain. 


8 Other Unresolved Issues 
8.1 X-ray Flaring 


The RXTEi PCA count rates (Fig. 2) show occasional brief maxima, which become 
dramatic in the year before the X-ray minimum . Farly in 1 997 1 shibashi noticed that 
the brightest X-ray ''flares" tended to recur every 85 days [50] . This periodicity was 
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Fig. 5 Possible orbit of companion star. The orbit has an eccentricity e = 0.9 with semi-major axis 
a 15 AU. Dots along the orbit show CHANDRA (CXO) observations, while images show the 
projected density in the orbital plane of the interacting winds of ij Car (high density; yellow-red) 
and the companion (lower density, blue-black) in one type of 3-d model (701. The curved arrowed 
line shows the best guess of the orientation of the line of sight projected onto the orbital plane. 


statistically doubtful, but it correctly predicted the next two major flares in July and 
September 1997. and perhaps another in March 1998 after the X-ray minimum [511. 

X-ray spectra show that they represent increased values of the volume emission 
measure [9, 49, 8], though some hardening of the spectrum during flares has been 
noted [ 14, 40|. The evidence so far favors a change in the hot colliding-wind plasma 
as the cause, rather than sporadic decreases in intervening absorption [9, 49). The 
physical mechanism producing the X-ray flaring is not understood. In a simple qual- 
itative model, localized density enhancements (clumps, or structured entities such 
as co-rotating streams [66 [) in the wind eventually encounter the wind-wind shock 
and provide a brief increase in the amount of shocked X-ray emitting material. The 
flow timescale through the shock is ffl 0W ~ a/v n ow where a is a characteristic size 
(thickness?) of the shocked region and pg ow is presumably a few hundred km s 1 
based on the wind speeds, or roughly 0.1 to 0.4 A1.J per day. The observed 10-30 
day flare durations are consistent with emission-region sizes of several AU [ 65 1 . 

The nature of the 85-day “quasi-period” is still not understood. A recurrence 
timescale of 85 days seems reasonable for either rotation or pulsation of the primary 
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star; and in some models it might vary from cycle to cycle if each spectroscopic 
event alters the star’s outermost structure [31 ]. But the interval between flare peaks 
varies with X-ray cycle phase, decreasing to as little as 20 days just before X-ray 
minimum, qualitatively consistent with expectations based on flow travel times |24|. 




Fig. 6 Timing of major and minor X-ray flares in 1997-1998 according to Ishibashi [49|. 


82 Cycle-to-Cycle Variations 

X-ray monitoring also revealed significant variations in the average X-ray flux level 
from cycle to cycle. Figures 2 and 7 shows RXTE count rates in the three 5.5-vear 
cycles.The first and third cycles had very- similar flux levels . but cycle 2 was signifi- 
cantly brighter. In a standard colliding-wind model the flux depends solely on orbital 
parameters, stellar mass-loss rates, and wind speeds (Section 3). Among these the 
mass-loss rate is probably most likely to change without perceptibly affecting the 
shape of the X-ray spectrum. 

Cycle-to-cycle changes have also been observed in other wavelength bands. He I 
10830 A emission showed a secular decrease for several cycles [16], centimeter- 
wavelength radio fluxes have been brighter in 2004-2008 than they were in 1998- 
2003 (S. White, priv. comm,), and the profiles of the extremely bright Ha and H/l 
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emission lines conspicuously differed between the 1997 and 2003 events |32|. Alto- 
gether, these developments likely indicate some undiagnosed change in the primary 
wind, although variations in other wavebands may be more complicated to inter- 
pret because they can be heavily influenced by secular changes in the Homunculus 
produced by expansion of the nebula and/or the creation or destruction of dust. 

The most recent X-ray minimum which began on 16 January 2009, Figure 7, 
was consistent with the phasings of the two previous minima. However, the X-ray 
recovery began on 18 February 2009, one month earlier than expected based on the 
duration of the two previous X-ray minima. During the 2009 recovery the X-ray 
hardness increased markedly, similar to the spectral change which occurred at the 
end of the previous two minima, although the maximum hardness observed in 2009 
was less than the maximum observed after the previous recoveries. This implies 
some fundamental change in the properties of the colliding wind emission. This 
may involve a decline in q Car’s mass loss rate [55,63, 15,64| 2 . though changes in 
the wind properties of the companion cannot be completely ruled out at this point. 


Phase 

- 0.05 ' 0.00 0.05 



Fig. 7 Comparison of three X-ray minima observed by RXTE. The start of the 2009 minimum 
(16 January 2009) was consistent with the ephemeris derived from the previous two minima. The 
bottom axis shows the Julian Day Number corresponding to the 2009 minimum. The X-ray flux 
recovery occurred about 30 days early compared to recoveries from the previous two minima. From 
Corcoran et at. [ 15]. 


2 See Figure 10 in the Humphreys and Martin chapter. 
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The X-ray minimum may be caused by an increase in the intervening absorption, a 
decrease of the intrinsic emission, or both; in either case the maximum temperature 
is observed to be fairly constant |52, 40] . All models so far agree that the event 
happens near periastron and, if an eclipse plays a leading role, the eclipsing agent is 
the primary wind and not the primary star itself. 

But several pieces of the puzzle do not easily fit into a simple eclipse model. 
The duration of the X-ray minimum is much longer than predicted in simple eclipse 
models, though more realistic wind models do better jobs in matching the eclipse 
duration [74, 70|. If the X-ray minimum really occurs near periastron, simple mod- 
els suggest that the intrinsic X-ray luminosity should be in excess of 10 35 ergs s' 1 
at E > 5 keV at that time since L x ~ 1 jd. As noted above, hiding high-energy emis- 
sion this bright requires an enormous wind column density which is about a factor 
of 100 higher than the observed column density. Another problem is that, during the 
minimum, a fully hidden X-ray source should produce a fluorescence iron line w ith 
an equivalent width about a factor of 6 larger than the observed value [40). Also, the 
X-ray minimum shows two distinct phases (a deep stage and a shallow stage)[40|, 
which is not expected in simple eclipse models. 

Several non-eclipse ideas have been proposed. A spectroscopic event may be 
the sporadic ejection of extra material by the primary star, perhaps triggered by the 
companion star’s close approach [100, 26, 11,61], This causes various effects, in- 
cluding a sudden rise of X-ray absorption. The latitudinal structure of the primary 
wind may suddenly change, increasing the wind density at equatorial latitudes where 
the secondary star is probably located [81]. Near periastron, the wind-wind shocks 
may disintegrate, collapse, or at least become very complex because of instabilities 
that arise when radiative cooling dominates [30,61 [.This would destroy most of the 
highest-temperature regions, thereby replacing the 2-10 keV emission by unobserv- 
able soft X-rays. Breaking the shock apart by accretion of part of the primary wind 
onto the companion star [83] could also reduce the colliding wind emission near pe- 
riastron. Near periastron, the secondary star might move into a dense structure such 
as a disk around the primary. These possibilities are not all mutually exclusive; in 
fact the first, second, or fifth of them would make the third more likely. Nor are they 
inconsistent with an eclipse, which, based on the likely orbit orientation, is likely to 
begin shortly thereafter [61 [. 

The near-constant slope of the continuum above 5 keV is evidence that the high- 
est temperatures are not reduced during periastron passage [40]. Finally, the strange 
He II A 4686 emission, which peaked in a dramatic way during the 2003.5 event, 
seems directly pertinent to the X-rays for two reasons: it rose quickly in brightness 
during exactly the lime when the 2-1 0 keV brightness was fading, and most likely it 
requires a remarkably increased flux of soft X-rays below 0.5 keV, which of course 
cannot be observed directly [61]. 
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9 Future Prospects 

X-ray observations provide unique direct information on the mass flows in the t; Car 
system and indirect information about the system and stellar parameters. The recent 
X-ray minimum in January 2009 has been observed by a fleet of X-ray observatories 
including KXTE, CHANDRA, XMM, Swift and Suzaku with a primary goal of con- 
straining the flow dynamics and determining the nature of the X-ray minimum. In 
the coming years high resolution X-ray spectra with X-ray calorimeters on observa- 
tories like the Astro-H satellite [90] should resolve the He-like iron K-line triplet at 
6.7 keV for the first time and provide well-defined profiles of the the H-like iron line 
near 7 keV. These observations will be important since the high formation temper- 
ature of these lines mean that they must be formed near the stagnation point of the 
wind-wind collision; tlius the line centroids can in principle be used to trace the or- 
bit of the stagnation point, which can in turn be directly related to the radial velocity 
curve of the companion. This would enable us to dynamically constrain the system 
mass function and provide key answers about the evolution of the binary. Resolv- 
ing the He-like iron K line triplet also provides important density and temperature 
diagnostics of the shocked gas near the stagnation point. 

Even higher-energy images of hard X-ray and Gamma-ray emission from t) Car 
have recently been reported[59, 89, 2], using modem observatories like the Interna- 
tional Gamma-Ray Astrophysics Laboratory (INTEGRAL), the Italian-led Gamma- 
ray satellite observatory Astro-rivelatore Gamma a Inmagi ni LEgge.ro (Agile), and 
the Fermi Gamma-Ray Space Telescope . As history shows, high-energy high-spatial 
resolution observations are extremely important, and revolutionary technology like 
the proposed space-based X-ray interferometer] 3 8] which could reach resolutions 
of micro-arcseconds at X-ray energies, holds out the hope of directly measuring the 
geometry of the colliding wind shock as it moves around the orbit. 

These observations will provide much of the the crucial information (system 
masses, stellar parameters) needed to fully understand the evolution of t| Car. Look- 
ing even further down the road. X-ray observations will play an important role in 
understanding the eventual core-collapse of the star and the formation of a black 
hole and its associated supernova remnant. 
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